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Abstract. - Using the numerical unrestricted Hartree-Fock approach, we study the ground state of a two- 
orbital model describing newly discovered FeAs-based superconductors. We observe the competition of a 
(0, n) mode spin-density wave and the superconductivity as the doping concentration changes. There might 
be a small region in the electron-doping side where the magnetism and superconductivity coexist. The super- 
conducting pairing is found to be spin singlet, orbital even, and mixed s^y + d^2 ,,2 wave (even parity). 



The newly discovered FeAs-based superconductors [1,2] 
have attracted lots of experimental [3-9] and theoretical [10- 
39] interests. Experimentally, the superconducting transition 
temperature, up to now, can be as high as 56K [4-6]. Spin- 
density wave order of (0, n) mode was observed in the parent 
compound LaOFeAs, but vanishes at high temperature (above 
150K) and large doping region [4,9]. However, unlike the 
cuprate high-Tc superconductors whose parent compound is 
an insulator, LaOFeAs is a semimetal [4, 10]. The observed 
magnetic dependence of the specific heat as well as the results 
of nuclear magnetic resonate suggest the presence of gapless 
nodal lines on the Fermi surface [3]. Theoretically, the transi- 
tion temperature estimated based on the electron-phonon cou- 
pling seems unlikely to explain the observed superconductiv- 
ity, thus suggesting these materials might be unconventional 
and non-electron-phonon mediated [16]. The local-density- 
approximation (LDA) calculations show that the density of 
state near the Fermi surfaces of the parent compound LaOFeAs 
are dominated by iron's 3d electrons [10-15]. These obser- 
vations imply that the multi-orbital effects play a key role in 
these new family of high-Tc superconductors. Though irons in 
LaOFeAs have five orbitals and ten bands, some groups pro- 
posed that it is sufficient to consider only a few of them, say d^-r 
and dy, orbitals for a minimal two-band model or may include 
dxy orbital for a three-band model, to reproduce qualitatively 
the LDA Fermi surface topology [20-29]. The spin-density- 
wave (SDW) of (0, n) mode of the parent compound LaOFeAs 
has been interpreted due to the superexchange interaction be- 
tween the next-nearest neighboring sites, which leads to an 



effective Ji - J2 model defined on a two-dimensional lattice 
[18, 25, 30, 31, 35]. Alternatively, the SDW mode was also at- 
tributed to the nesting of Fermi surface [4, 12,21,27,28,33,34]. 
On the other hand, despite of the pairing mechanism being un- 
clear, some groups [36-39] have addressed the classification of 
various superconducting states of the systems with two orbitals 
based on group theory. 

In this paper, we study the competition between supercon- 
ductivity and magnetism, as well as the superconducting pair- 
ing symmetry in the FeAs-based materials. Based on a recently 
proposed two-orbital model, we introduce the superconducting 
pairing into the Hamiltonian from which various pairing possi- 
bility can be constructed. Then we study the ground state prop- 
erty by using the unrestricted Hartree-Fock approach. We start 
with random values for all possible order parameters and let the 
self-consistent equations converge to the final solution. We find 
that for the parent compound, there exists a (0, n) mode spin- 
density-wave. When doping is introduced, either by electron 
or hole, superconductivity appears with the disappearance of 
the magnetism, and there is a small overlap in electron-doping 
region where the magnetism and superconductivity might coex- 
ist. Therefore, these is a competition between the two orders. In 
the superconducting region, a careful scrutiny reveals that the 
pairing symmetry of the superconductivity might involve spin 
singlet, orbital symmetric, and mixed St,,H-dj.2_^,2 wave (even 
parity). 

In order to simplify the calculation, we adopt the two-band 
model of the Fe subsystem defined on the principle plane, as 
suggested by Raghu etal [23], as a starting model for discussing 
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Fig. 1 : (color online) The schematic illustration for the hopping pa- 
rameters and the attractive interaction parameters of the two-orbital 
dxT. — dy^ model on a square lattice. 



the interplay between magnetism and superconductivity. While 
there are arguments as to whether the minimum model should 
consist of three orbitals, we expect no qualitative difference to 
appear within our calculation framework. There are two ''dxz, 
dy,'" orbitals on each site. The left plot of fig. 1 shows all pos- 
sible hoping processes between neighboring sites and orbitals 
based on the tight binding approximations. Phenomenologi- 
cally, we use the following mean-field Hamiltonian, 



H - Ha + Hsc + Hi, 

Hq ^ 2_j Kit [^*^0 + yi<'^^ + f*'^3] ^^ko 



(1) 

(2) 
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The physical meaning of all parameters and operators in the 
Hamiltonian are presented in order. Hq is the noninteracting 
term [23] with 



(CUk)X'jk)), 



(4) 



£k = -(f 1 + f2)(cos k^ + COS ky) - 4f3 COS k-c COS ky - fi{5) 
jk - -4f4sinfci-sin^y, (6) 

^k = -(f 1 - f2)(cos k^ - cos ky), (7) 

and a,yS - a(xz),b(yz) label orbitals, and r is the Pauli ma- 
trix for xz and yz orbitals [The superscript a(b) is specified to 
"dxzidyzY' orbital hereafter]. Hj is a Hartree-Fock decomposi- 
tion of the on-site interaction term, result from 



Hi = f/X<T«a 



i./3 
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Fig. 2: (color online) The superconducting order, the magnetism order, 
and doping as function of the chemical potential. 



in which Jh denotes the Hunds rule coupling, U{U') the on- 
site Coulomb interaction between electrons on the same (dis- 
tinct) bands, and 5^" = Cf!cf|,5f = (/. - r^n)l2,j^^ = 



rlf' = (/ - l)/2. Both S and 77 satisfy SU(2) Lie alge- 
bra, and U - U' + Jh [40] . Local order parameters nt: and p! 






are defined as 



4 - ^(<<T>-<<i>)' 



^ i.a' 



(9) 

(10) 



Respecting to the fact that underline mechanism for super- 
conducting properties observed in cuprates is controversial, we 
are not in the position to discuss what mechanism is for the 
FeAs-based superconductivity. Instead, our purpose is to ad- 
dress the interplay between superconductivity and magnetism 
so we will start from the assumption that electrons on neigh- 
boring sites pair in the following form 



•^^ / J 1,6,0- i+6,cr i,-(T I i,S,o 



i,6,p,a- 



(11) 



-V y Ut Ci„ C"^ + h.c.-\At 

/ i 1,6' .(7 t+d ,(T r-iT I 1,6 ,cr 

i,d',a,p,a' 

where V accounts for the attractive interaction between elec- 
trons on neighboring sites, S(6') denotes the (next) nearest 
neighboring sites [fig. 1 (Right)]. 

The mean-field Hamiltonian is quadratic and can be diag- 
onalized numerically by solving self-consistent equations fol- 
lowing standard procedure. In our numerical calculations, we 
set fi = -1, f2 = 1-3, f3 = f4 = -0.85, and Jh = 1.4 ~ 2.0 
in unit of |fi|, according the LDA calculations, and vary other 
parameters to explore the ground-state prosperities. Solving 
such multi-variable nonlinear equations is a very complicated 
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and time consuming task because of many quantities under it- 
eration: we have 8 independent A^^^^s, 16 A^'^, ^s, 4 A^ s, and 

4 density parameters p.^s, for each site, and for a L x L-site 
system, so the total number of iterated quantities is 321?. The 
reason for us to perform the unrestricted Hartree-Fock calcu- 
lations, rather than the restricted ones (where one would have 
much reduced number of self-consistent variables), is to obtain 
unbiased ground state pairing and magnetic ordering patterns. 
To address the competition between the magnetism and su- 
perconductivity. We first define the superconducting and mag- 
netic order parameters in ^-space as follows 



ji 
nil. — — ^ > e'-' m.. 



Va^^ 



(12) 



(13) 



Then we calculate the magnetic order and superconducting or- 
der parameter defined as. 



= -yki 



(14) 
(15) 



Fig. 2 shows the main results of the doping ratio, magnetic or- 
der, and superconducting order as a function of the chemical 
potential for a typical parameter set U — 2>A, V - 0.9, and 
Jh - 1-5. Note that paring driving force usually comes from 
a second order process so it is expected that V is smaller than 
any other renormalized Coulomb interactions. We firstly ob- 
serve that there is a flat region in the doping as the chemical 
potential varies. Though the position of the flat is a little above 
the half-filling for finite size lattices we studied, it will tend to 
the half-filling as the system size increases. From fig. 2(b), we 
notice that the magnetic order is nonzero around the half-filling. 
This observation is consistent with experimental results that the 
parent compound shows magnetic order without doping. The 
dominant magnetic order is (n, 0), or (0, n) SDW. Once elec- 
trons or holes are doped into the parent compound, the mag- 
netic order is gradually suppressed to zero, and superconduct- 
ing order appears, as shown in fig. 2(a). At half-filling, the 
superconductivity is absent. In experiments, the parent com- 
pound is a poor metal without superconductivity. Therefore, 
our numerical results are consistent with the experimental re- 
sults qualitatively. As electrons are doped into the material, the 
superconducting order appears at a certain doing concentration. 
As the doping concentration becomes larger, it will disappear 
again. Similar behavior also occurs in the hole doping region. 
One of the most interesting observation is that, in the electrons 
doping region, both the magnetism and superconducting orders 
do not vanish at a small region. Therefore, the magnetism and 
superconductivity are not exclusive with each other So though 
there is a competition between two orders, the magnetism and 
superconductivity might coexist in a very narrow region. 

Next we check the magnetic order and pair symmetry of the 
superconductivity. The order parameters defined in A:-space, i.e. 
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Fig. 3: (color online) (a) (0, tt) spin density wave at half-filling. 
(b)Pairing structures of the (next) nearest neighboring sites in electron 
doped region, (c) Pairing structures of the (next) nearest neighboring 
sites in hole doped region. Here A,, i = 1,2,3,4 denote simply the 
pairing amplitude which will be shown in main text. 



eq. (12) and eq. (13), can help us to learn the details of mag- 
netic order and pairing symmetry. In this work, we consider 
three typical cases below. 

Case A (undoped compound): We consider the magnetic 
property of the parent compound for the parameter set 



U = 3.5, V = 1.0, Jh = 1.5, ju = 3.3. 



(16) 



Here, we set the initial value of local magnetization at random 
and let self-consistent equations converge. We observe that the 
dominant magnetization patten in k space include either (0, n) 
or (tt, 0). Fig. 3(a) shows the magnetic ordering in both the real 
and A:-space. In the real space, the compound shows a stripe- 
like phase which represents the spin collinear order The local 
magnetization are ±0.227 for light and dark gray respectively. 
This observation is consistent with the experimental results of 
the magnetic properties in undoped case. If we make a Fourier 
transformation, we find a spin-density wave of (0, n) mode in 
k space. Therefore, irik - 0.227 at (0, ±7r), and zero elsewhere. 
Here, we would like to point out that (:7r, 0) mode can also be 
observed even we started iteration with initial conditions con- 
sist with of other order parameters. Moreover, our numerical 
simulations show that the magnetic order dominates in the low- 
doping region. 

Case B (electron doping): Electron doping typically occurs 
at 

f/ = 3.5,y= 1.0,7h = 1.5,yu = 3.7, (17) 

where the solution converges at 13.8% electron doping. In 
this large doping region, the anti-ferromagnetism of the par- 
ent compound are suppressed by superconductivity. The intra- 
orbital pairing structures of the (next) nearest neighboring sites 
are illustrated in fig. 3(b) schematically. In fig 3(b), Ai = 
<C-;c:l,,) - 0.04824, A2 = (S^Ct,,^ = 0.05176 and 
A3 = {Cf^C"l^^y^) = 0.03584. While the inter-orbital pair- 
ings are found to be zero within numerical accuracy. Hence 
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from numerical results, the nearest-neighbor intra-band pairing 
belongs to dj^2_^2-wave symmetry (cosk^ - cos^y), while the 
next-nearest-neighbor intra-band pairing looks like an extended 
s^y-wave (cos kj, cos ky), and there is no inter-band (next) near- 
est neighboring pairing. In addition, the superconducting order 
parameters in k space obey the following relation 



(cf c«l ) = (c''! cf ). 

\ k,(T -k,-(T/ \ -k,(T k-o-f ■ 



(18) 



which concludes that the pairing symmetry involves spin sin- 
glet, orbital symmetric, and mixed s^^ -i- dj.2_,,2 wave (even par- 
ity). 

Case C (hole doping): Hole doping typically occurs at 



f/ = 3.5,y= 1.0, 7^= 1.5,1/ = 2.8 



(19) 



where the solution converges at 11.5% hole doping and with- 
out spin density wave. The first two graphs of fig. 3(c) show 
that the (next) nearest intra-band pairing are still mixed s,xy 
+ d^2_j,2 wave. The pairing amplitude are Ai = 0.04525, 
A2 - 0.01864 and A3 = 0.03444. Moreover, unlike case B, the 
next nearest inter-band pairing interactions are not too small 
to ignore. As illustrated by the last two graphs in fig. 3(c), 
the signs of the order parameters in real space denote the d- 
wave inter-band pairing structures. The pairing amplitude is 
A4 = {C"\C'''^ , ) = 0.00185. We also notice that the relation 

{cflC^l^o) = {C%C''^!-c-) still holds in the hole doping re- 
gion. Therefore, the pairing symmetry are the same as in case 
B. 

In summary, applying the numerical unrestricted Hartree- 
Fock approach to a two-orbital model, we have studied the 
competition between magnetism and superconductivity, as well 
as superconducting pair symmetry in the newly discovered 
FeAs-based superconductors. We found that there does exist 
competitions between magnetic and superconducting orders. 
Despite of this, the magnetism and superconductivity may still 
coexist in a small region. In order to find the possible mag- 
netic order and pairing symmetry. We further studied three dif- 
ferent cases, corresponding to undoped, electron doped, and 
hole doped, respectively. We found that, around updoped re- 
gion, the parent compound shows a spin-density wave of (0, n) 
mode, while in both the electron and hole doped region, the su- 
perconducting pairing symmetry involves spin singlet, orbital 
symmetric, and mixed s^^ + dj^i^yi wave(even parity). We no- 
ticed also that the result is somehow different from that of nu- 
merical renormalization calculation [33]. This might be due to 
the different range of interactions, as well as different approach 
used. Finally, we remark that although we used the two-orbital 
model, suggested by previous studies, as the starting model in 
our calculation, conclusions drawn should not be altered quali- 
tatively when more orbitals are included in our approach. 
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